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Abstract. A family of random matrix ensembles interpolating between the 
GUE and the Ginibre ensemble of n x n matrices with iid centered complex 
Gaussian entries is considered. The asymptotic spectral distribution in these 
models is uniform in an ellipse in the complex plane, which collapses to an 
interval of the real line as the degree of non-Hermiticity diminishes. Scaling 
limit theorems are proven for the eigenvalue point process at the rightmost edge 
of the spectrum, and it is shown that a non-trivial transition occurs between 
Poisson and Airy point process statistics when the ratio of the axes of the 
supporting ellipse is of order n"^/''. In this regime, the family of limiting 
probability distributions of the maximum of the real parts of the eigenvalues 
interpolates between the Gumbel and Tracy- Widom distributions. 

1. Introduction of the model 

1.1. Introduction. The local statistics at the edge of the spectrum are of funda- 
mental interest in the study of random matrix models. For the Gaussian unitary 
ensemble (GUE), defined as the probability measure 

(1-1) - exp{-^Tri7^}di/ 

2 

on the space "Hn ^ K" of Hermitian matrices, where AH denotes Lebesgue mea- 
sure, the celebrated work of Tracy and Widom, [14], shows that the distribution 
of the largest eigenvalue, properly rescaled, converges to Ftw, the Tracy-Widom 
distribution, which is expressible in terms of a solution to the Painleve II differential 
equation. Their results have been generalized to other models and it is conjectured 
that Ftw is the universal limiting extreme value law for large classes of random 
matrix ensembles. This universality conjecture is supported by the fact that the 
Tracy-Widom distribution has appeared in a number of seemingly unrelated con- 
texts, see for instance [1], [8] and |9]; in some sense it is a natural extreme value 
distribution. 

In contrast, the classical extreme value theorem for a set of independent identically 
distributed (iid) random variables, say Xi e A^(0, 1/2), states that the distribution 
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of their appropriately shifted and rescaled maximum converges to the Gumbel dis- 
tribution, Fc{x) := e~'^ , 

ma.Xi<i<n{Xi} — Cn 



lim 

n — >c< 

where 
and 



< t 



log(47r logn) 
4 Vlog n 



2 Vlog n 

Although the edge behaviours of these models have received much less attention 
in the literature, there are random matrix ensembles without symmetry conditions 
imposed so that the Hmiting spectral density is no longer supported on a line, but 
some other set in the complex plane. The simplest example of such a model is the 
classical Ginibre ensemble of matrices with independent centered complex Gaussian 
entries of variance n~^, introduced in [6], for which the asymptotic spectral distri- 
bution is uniform in (and vanishes outside) the unit disc. There are several possible 
notions of "extremeness" that might be considered in this model; Rider [12j shows 
that the spectral radius, that is, the maximum absolute value of the eigenvalues, is 
asymptotically Gumbel-distributed after a shifting and rescaling much the same as 
for the case of extremes of independent random variables. Intuitively, this can be 
understood as reflecting the fact that the few eigenvalues with largest modulus do 
not in general lie close to or feel the repulsion of each other, and thus behave as if 
they were independent. Rider's method rehes heavily on the radial symmetry of the 
problem; a closely related problem which generalizes easier to the family of models 
to be considered here is that of determining the distribution of the maximum of the 
real parts of the eigenvalues instead of their absolute values. 

1.2. The ellipse ensemble. The fundamental difference between the edge statistics 
of the QUE and those of the Ginibre ensemble corresponds to the distinction between 
random matrix and Poisson statistics which has attracted a great deal of attention 
from the mathematical physics community in connection with quantum chaos. In 
particular, it is illuminating to study families of models that exhibit a transition 
between the two types of behaviour, see e.g. [lOj. Based on the observation that 
a Ginibre matrix A can equivalently be considered as a sum A — {Hi + iH2)/V^ 
where Hi and H2 are two independent GUEs, a natural family of random matrix 
ensembles interpolating between the GUE and Ginibre ensembles can be defined as 
the distribution of 



^.y(i±i)„, + .y(i_i)„,. 

for T e [0, 1]. The parameter t can thus be thought of as determining the extent 
to which A fails to be Hermitian. In j^, the authors define this ellipse ensemble 
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explicitly as the probability measure on the space Mn of complex nxn matrices 
given, for T < 1, by 

(1.3) dPliA) ^ (^-—^^==^ exp|-^^^Tr(AAt-rReA2)|dA 

where dA is Lebesgue measure on A^„ ^ R^"^. It is easy to see that the cases t — 
and T — > 1 correspond to the Ginibre and Gaussian unitary ensembles respectively. 
The successful analysis of this model depends crucially on the fact that, for every 
T e (0, 1), the induced eigenvalue measure 

(1.4) dp;'(Ci,...,C„) 

=^ n ia--opexp|--^x:(ioi^--i^«ci)|d%, 

^ l<]<k<n y ^ ' j=l J 

where ZJ^ is a normalizing constant, can be considered as the distribution of a 
point process Z^^ on (= C) with determinantal correlation functions given by the 
correlation kernel 

(1.5) x;((ei,m),(6,??2)) 

Here Cj = + iVj, and hk are the normalized Hermite polynomials, that is, the 
orthonormal polynomials with respect to the measure dx on R. (See Section [2?T] 
for definitions related to determinantal point processes.) Furthermore, the rescaled 
Hermite polynomials appearing in l|1.5p are orthogonal in the whole complex plane 
with respect to the corresponding exponential weight. The orthogonal polynomial 
techniques employed in the analysis of the GUE therefore become available for the 
ellipse ensemble as well, although considerable technical difficulties remain, since 
the Christoffel-Darboux formula no longer appHes to simplify the sum of products of 
rescaled Hermite polynomials. It is proven in |4j e.g. that for fixed r ^ 1, the Hmiting 
eigenvalue density is constant in the ellipse {(x, y) : j (1 + r)^ + j/-^/ (1 — r)^ < 1}, 
where z — x+iy. The authors also discover a regime of "weak non-Hermiticity", when 
T = Tn = 'V — o? I {2n) for a new parameter < a < oo, in which a transition occurs in 
the local bulk statistics between the sine kernel behaviour of the GUE and the scaled 
correlation kernel of the Ginibre ensemble. In [5], it is shown heuristically that in 
this regime the local statistics near the rightmost edge of the spectrum are given by 
a correlation kernel which factorizes as the Airy kernel in the ^-variables multiplied 
by a Gaussian weight in the ry-direction; thus the model essentially behaves like the 
GUE at the edge in this regime. 

However, as the present work will show, there is another transitional regime, for 
(1 — r„) ~ n^^/^, when the edge statistics interpolate between those of the GUE and 
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the Ginibre ensemble. Let (t„ = n^^^y/{l — t„). It is proven in Theorem 12.51 that if 
cr„ — > 0, the rescaled eigenvalues near the rightmost edge of the spectrum converge 
to a point process on which can be interpreted as the Airy point process in the 
x-direction with each particle subject to an iid centered Gaussian displacement in 
the y-direction. Furthermore, the distribution of the maximum of the real parts 
of the eigenvalues converges to the Tracy- Widom distribution. The process is thus 
essentially one-dimensional and coincides with the corresponding point process at the 
edge of the GUE. (This is what is found in [5] for the particular choice (T„ oc If 
(T„ oo the edge point process converges to a Poisson process on with intensity 
7r^^/^e~^~^ , and in addition the maximum of the real parts of the eigenvalues 
converges in distribution to a Gumbel random variable. In the intermediate regime 
when an cr G K."*" a new non-trivial interpolating point process on arises in the 
n ^ CJO limit. It is again a determinantal point process with a rightmost particle 
almost surely, but the correlation kernel no longer factorizes. 

The edge behaviour in the various regimes can be interpreted as follows: When 
(T„ ^ oo, the supporting ellipse collapses to an interval of the real line so slowly that 
the imaginary parts of the eigenvalues near the edge are of order much greater than 
the spacing of their real parts. Therefore the eigenvalues are not close and do not 
interact, so in the appropriate scaling Hmit they behave as a Poisson process. On the 
other hand, when ct„ 0, the imaginary parts of the eigenvalues near the edge are 
negligible compared to the spacing of their real parts; hence they behave essentially 
like the particles in the Airy point process with independent fluctuations in the y- 
direction. In contrast, the case (t„ — > cr corresponds precisely to the imaginary parts 
of the eigenvalues being of the same order of magnitude as the spacings of their real 
parts, and there is a non-trivial interaction between eigenvalues even though their 
displacement away from the real axis is no longer negligible. 

The proof relies in an essential way on flnding, in the absence of a Christoffel- 
Darboux formula, a manageable representation of the sum of products of Hermite 
polynomials in (|1.5p . This is provided by the double integral formula of Lemma SHJ 
The asymptotics of the correlation kernel can then be analyzed by means of saddle 
point techniques. 

Although a very speciflc family of models, particularly accessible to analysis, is 
studied here, it is reasonable to conjecture, in analogy with the universality con- 
jectures for random matrix ensembles with purely real spectrum, that the type of 
transition and scale at which it occurs has a universal character that should be found 
in much greater generality for random matrix models in which the support of the 
spectral density collapses to a one-dimensional subset of the complex plane. Exam- 
ples of other non-Hermitian ensembles are provided by the model introduced in [15] 
and analyzed rigourously in [?]• 
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2. Formulation of results 

2.1. Preliminaries on determinantal point processes. For a comprehensive ac- 
count of the theory of point processes, see [2]. The material relating to determinantal 
processes can be found e.g. in [11]. 

Let A be a complete separable metric space. A Borel measure ^ on A is a houndedly 
finite simple counting measure if it takes non-negative integer values on bounded sets 
and ^({a;}) G {0, 1} for every a; G A. The set of all boundedly finite simple 

counting measures on A can be identified with the family of sets x = {xi\ C A such 
that Xi 7^ Xj \{i ^ j and n i?| < oo for every bounded B C A. Define the smallest 
cr-algebra on A/'(A) such that the mappings A/'(A) 9 /i i— > fi{A) are measurable for 
each Borel set A G A, making Af{A) a measurable space. A point process X — {xj} 
on A is a random element of Af{A) . 

Let A be a reference measure on A (in our case, A will be K or and A Lebesgue 
measure). If AT is a point process on A and, for some n > 1, there exists a measurable 
function p„ : A" R such that for every bounded measurable function (p : A" ^ M., 



(2.1) 



E 



,Xk„ 



^(6,---,CnK(Cl,---,en)d"A(C), 



then pn is called an n-point correlation function of X. (The sum on the left hand 
side of (|2.ip is over all n-tuples of distinct points in X, including permutations.) It 
can be shown that if X has correlation functions p„ for every ti > 1 (and the product 
over an empty index set is 1 by definition), then 



(2.2) 



E 



l[{l + cj){x,)) 



n=0 j=l 



for every bounded measurable function (f) with bounded support. (Here the product 
on the left hand side is over all particles in the point process; since there are only 
finitely many particles in each bounded set the product is finite for each realization 
of X). 

A point process for which all correlation functions exist and are of the form 

Pn(6, • ■ • ,60 = det {K{^„Q)l^^^ 

for some measurable function K : A^ ^ C {a correlation kernel) is called a determi- 
nantal (point) process. The correlation kernel K is said to be Hermitian if 



if(6,ei) = ^(a,6) 

for all (a,6) e A2. 

Weak convergence in N[A) of a sequence {Xn} of point processes to a Hmit X is 
equivalent to weak convergence of the finite dimensional distributions, meaning that 
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for any finite family {Aij^L^ of bounded disjoint Borel sets in A and integers rrii, 

Jim^P„ [\Xn nA\^ m,]j = P [\X nA\^ m,]j . 

Consider a determinantal point process Z = {zj} — {{xj,yj)} on M^. If there is 
a i e M such that E [\Z n ((t, co) x M) |] < co, Z is said to have a rightmost or last 
particle almost surely. The last particle distribution function, F, of Z is then defined 

as 

(2.3) F(t) =P[|Zn ((i,oo) X M)| = 0] . 

To prove the limit theorems of this paper we will appeal to the following lemma, 
which is a fairly standard result. For completeness, a proof is given in an appendix. 

Lemma 2.1. For every positive integer n, let Zn be a determinantal point process 
on K.^ with Hermitian correlation kernel if„. Suppose that Kn{C,i, C2) (2) as 

n ^ 00 for every (Ci, (^2) € and that there is a measurable function S : ^ M 
with -B(C) > |Jsr„(^,^)| for every n, such that for each ^0 G K 



/ i3(C)dC - < 00. 

"'(fn.oolxR 



'(?o,oo)xl 

Then Zn converges weakly as n 00 to a determinantal point process Z with corre- 
lation kernel K , Z„ has a last particle almost surely with distribution 

(2.4) Fn{t) = £ / det(if„(0, a))J,fc.id'-C, 

r=0 "'((t,oo)xR)'' 

and Fn converges to the last particle distribution F of Z as n ^ co. 

2.2. The interpolating process. The Airy kernel is the function Ka : R 
defined by 

^ Ai(gi)Ai'(6)-Ai'(ei)Ai(6) 

J<A[tl,t,2) = T T 

if $1 ^ $2 and by continuity if ^1 = ^2- Here Ai : M ^ M is the C°° Airy function, 

Ai{t) = — I e3"'+*"*dw e M, 

and 7 is the contour i t + i(5 for some 5 > Q, independently of the choice of 5 by 
Cauchy's theorem. We will use the double integral representation 

of the Airy kernel. It is well known (see [lij and [3]) that Ka defines a trace class 
integral operator on L^(t, 00) for any real t and that it is the correlation kernel of a 
determinantal point process Xa on M, the Airy point process. Xa has a last particle 
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almost surely, the distribution of which is known as the Tracy- Widom distribution, 
given by the Fredholm determinant 

Frwit) := det(/ - KA)L-(t,oo) = / det(AAte, ^,))t=id'e 



k=0 



Alternatively, 

(2.6) FtvfW = expj-^ {x-t)q'^{x)dx^ , 

where q is the unique solution to the Painleve II equation 

q"{t)=tq{t) + 2q{tf 

with the asymptotic behaviour 

q{t) Ai{t) as f — > oo. 

We now construct a point process Za on which, intuitively, to each realization 
{xi} of Xa associates a realization {{xi, yi)} of Za, where the yi are picked indepen- 
dently from A''(0, 1/2). More precisely, for every fc > 1, Za should have correlation 
functions 

fc 

Pfe(Ci,---,Cfc) 



:= det(i^Ate,6))t=in-^ 

i=l 

7^ 



(2.7) = det I ^ if^(ei, ij) 

) i,j=i 

making Za a dctcrminantal process as well. 

A Poisson process Zp on with intensity 7r~^/^e~^~'' can be viewed as a de- 
generate determinantal point process with correlation kernel 

(2.8) Mp(Ci,C2)=<5c,c,^ 



7r-V2e-^i-'?? if Ci =C2 
\ otherwise. 

To stress the analogy with the construction above, define 
(2.9) Mpi(Ci,C2) = y= i^p(a,6), 

where Kp{^i,^2) = ^iii2^^^^ is the correlation kernel of a Poisson process Xp on 
M with intensity e~^. Xp too has a last particle almost surely, which is Gumbel- 
distributed, since 

det{Kp{^„^,))l^^,d-^ 
E ^ / n e-^M"^ = E (e-r = i.,(,). 

n=0 •^(t.oo)"i=l „=0 '^^ 
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For any n > 1, 

det(Mp(C„a))7,fc=i =det(Mpi(0,a))"fc=i a.e. in R^", 
so Mpi can equivalently be chosen as a correlation kernel of Zp . Similarly, 

(2.10) Mp2(Ci,C2) = ^.nm T 

is also an equivalent correlation kernel of Zp. The following proposition sums up 
these observations. 

Proposition 2.2. The determinantal point processes Zp and Za on with corre- 
lation kernels Mp (or equivalently Mpi or Mp2) and 

(2.11) ma(Ci,c2) = - — j= — i^A(a,6) 

respectively, both have last particles almost surely, with distribution functions Fq 
and Ftw respectively. 

Proof. If iiT is a Hermitian correlation kernel of a determinantal point process on E 
satisfying 

/oo 
if(C,C)d^<oo, 

it is immediate that M((i,<^2) — 7r^^/^e"3(''i+''2)^(^^^ ^2) satisfies the correspond- 
ing conditions of Lemma IA.2I Since the correlation functions factorize by (|2.7p , 
the ?7i-variables may be integrated out in the expression l|A.2p for the last particle 
distribution. □ 

The following theorem asserts the existence of a family of determinantal point 
processes {Za}a>Q on R^ which, appropriately rescaled, interpolates between Za 
and Zp. 

Theorem 2.3. (Interpolating process) For each a G [0, 00) there exists a determi- 
nantal point process 

on with correlation kernel 

2 f r ^(fi'— '?2)^+f iJ''+iC2f — ^(<™+'7i)^+3«^+«4i" 
(2-12) d.d., 

where Q — (Cj,'7j) and is the contour t t-^ J jit) — t + iSj, independently of the 
choice of Sj > 0. 

Define the rescaled point process 



Za 

where 



'3 '"O' 

da ' be 

a 

V 6 log a ' 
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(6l0gCT)l/4^ 



3 log ^ ~ ^ log(6 log (t) — log(27r) 

The processes appropriately reseated, interpolate between Za and Zp in the sense 
that Zq = Za and Z^, converges weakly in A/'(M^) to Zp as a ^ oo. For any fixed 
a, Z„ has a last particle almost surely, with distribution function 



F,{t)=¥, [|Z,n((i,oo) xM)| 
= det(/ - Mcr)L-^((t,oo)xM) 



0] 



(2.13) 



r=0 



((i,co)xI 



det(Af.(C„a))-.=id^ed'-77, 



and furthermore 
(2.14) 



Fa{ca + a^t) Fa{t) as a ^ oo. 



Remark 2.4. To the author's knowledge, the processes Z^ have not been studied 
previously. In particular, we are not aware of any generalization to cr > of the 
representation (|2.6p . for the case cr = 0, of the Fredholm determinant l|2.13p . 

The proof is deferred to Section [31 

2.3. Edge scaling limits. This section contains the statement of the main result, 
giving scaling limits around the rightmost particle of the eigenvalue point process 
Z^ = {zj}"^i = {(a:^ , of the ellipse ensemble, defined formally as the 

measurable function 

mapping A G Ain to the set {zj}"^i of eigenvalues of A (with probability one, there 
are n distinct (complex) eigenvalues and we assume this to be the case) . Recall that 
ZI," is determinantal with correlation kernel KJ^" . 



Theorem 2.5. Let {t„}^]^ C [0, 1) be given and put an = v}!^ \J (1 — t„). For the 
choices of scaling parameters a„, 6„ and c„ specified below, define the rescaled edge 
eigenvalue point process 



zi- 



n/ J j = l 



and let 



F„-"(i)=P;" [^,I"n((i,(^) xK) 
maxi<j<„{a;j} - c, 



= 
< t 



be the last particle distribution of Z^' 
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(i) Suppose (Jn oo as n oo. Choose 

(2.15) a„ = fy2>'^~'^' 



V6 log cr„ ' 

^5/2 2/3 

^^ ic\ I --1/4 

(2.16) &n = T„ / 



(61oga„)i/4' 



anfl 



(2.17) dn = 2r„ + 5„ (^3 log (7„ - ^ log(6 log a„) - log(2^f3/4)^ , 

where f„ := (1 + t„)/2. T/ien Z,^" converges weakly to Zp and F^"-{t) 
converges to FQ{t) as n ^ oo. 
(ii) Suppose (T„ ^ ct e [0,cx)) as n oo. Choose a„ — ■nr'^l'^ , 6„ = cr„n^2/^ 
anrf c„ = (1 + Tn). T/ien Z^" converges weakly to Z^r and F^"-{t) converges 
to Fa{t) as n ^ oo. 

In particular, if a = 0, Z^ converges weakly to Za and F^" (t) converges 
to Frwit). 

Remark 2.6. Although the case t„ = has not been explicitly included in the 
calculations, a similar (simpler) argument can be made in this case (see the comment 
following l|4.10p ): as can be expected from the fact that the ellipse ensemble for 
constant r G (0, 1) is just a trivial rescaling of the Ginibre ensemble (see [4]), formally 
inserting r„ = in (|2.15p through (|2.17p indeed gives the correct limit theorem for 
the pure Ginibre ensemble as well. 

Remark 2.7. In view of the scaling Hmits in Theorems 12.31 and 12.51 (ii). the scaHng 
limit in Theorem l2.5l (i) is the only reasonable candidate (as long as t„ ^ 1); together 
these theorems in essence assert that Zn^"'' converges to Zp regardless of whether a 
tends to infinity with n or the n limit is taken first. 

3. Proof of Theorem 12^ 



The asymptotic notation defined in the first paragraph of Section [4] is used in this 
section as well, with n replaced by a. 



Proof of Theorem \2.3[ The existence of Zcr, and the fact that it has a last particle 
with distribution l|2.13p , follows by Lemma 12.11 from the convergence of Z^" when 
(t„)5^2 is a sequence such that 71^/^^/(1 — t„) ^ a as n ^ oo. This is established 
in the proof of Theorem 12.51 Note that Zq coincides with Za, by (|2.5p . 

To prove weak convergence of to Zp using Lemma [2rT| we first prove point-wise 
convergence of the correlation kernels. From the definition of correlation functions 
(|2.ip . it is clear that the rescaled point process Za- is determinantal with correlation 
kernel 

(3.1) Ma{Ci,C2) aabaMaiica + aa^i,bar]i),{ca + aa^2,bar]2)). 
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We fix (Ci, C2) = ((^1) '7i); ($2, '72)) G and calculate the cr ^ 00 limit of Mo- by a 
saddle point argument. Let 



(3.2) 

and 

(3.3) 



1 



Uiu) = --{au + ri[y + + iCiU, 



9<t{v) 



where = + a,a^j and rf^ = ba-rij . Choose the saddle points uq of and vq of 
satisfying 



and 



Uo 



Vo = 



a cr* 



77^ , , 2irl^ 2C2V2 



O 



respectively. Note that 



(75/2 (log C7)3/4 

1 

a5/2(loga)3/4 



,>o)=-.^ + 2..o = -.^(l + C9(-37^ 



)l/4 

and similarly f'^iuo) = —a^ + 2iuo. After the change of variables s = a{u — Uo), 
t = a{v — Vq), and choosing Si = Im(uo) and 62 = Im(wo) so that the new contours 
of integration become the real axis, the rescaled kernel can be written 

of a (s/o-+Mo )+5<T (t/a+VQ ) 



(3.4) 



M<.(Ci,C2) = 



aab. 



47r5/2 



o-^ yiM2 i{t/a ■ 



-dsdt. 



■ + uo + s/a + Vo) 

The main contribution comes from near the saddle points; let / = (— ro,ro) for some 
K ro « cr^/^ (log cr)^/^ and note that 



(3.5) 



IL 



o/a (s/o'+«o)+Sa- {t/(7+Vo ) 



47r5/2(T2 JJmx(r\i) i{t/a + Uo + s/a + vq) 



dsdt 



< 



aaha 



g/,(s/<7+«o)ds 



/ 



Now 

ga{t/(J + V(i) 

so putting 



= 9a{vo) + t^ 



1 

'2 a 



Wo 
2 



/ 7 *'f^n 

t -avo + V2 + — 



a 



3a3' 



CT3/2(log (7)1/4 
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and 



gives 



£2 = Re I -cTVo + + — + =0 



(loga)i/4 



t3/2 



g3<,(t/cr+t)o)^^ 



< 



,9<t(i'o) 



By the same argument for fa, and l|3.5p . it follows that 



(3.6) 



ofir{s/rT+Uo)+g^{t/lT+Vo) 



-dsdt 



47^5/2^.2 J J^2\j2 i{t/a + uq + s/(t + vq) 

<Ci- ' ' 



Since Im(cri + ct^uq + + (t'^vq) x 2c„ > 0, we can write 

1 ,ip(a(t+s)+(j2(„o+.uo)) 



< c- 



ro 



Therefore, using p.6p and Fubini's theorem, i|3.4p becomes 



dp. 



M.(Cl,C2) = 



aabo 



p/ct (s/ct + UO ) +ff<, (t/o-+«0 ) 



47r5/20-2 77^2 i(t/tT + Uq + s/a + Vq) 



dsdi + o(l) 



4^5/2 



oo 

e 

-aabaC^"'-^"''^^''-""^ f°° 



(3.7) 



2^3/2 

—a 5 gg^('"o)+f^{uo) 



where the last equality follows from the estimate 

1 



+ 0(1), 







(3.8) 



Pa + ha 



< c 



1 



obtained by integration by parts; in this case (in/ -Ja^ i< \J 6 log a ~> 1. Factors of 
the form exp{F(fi, 77^) — F(^2: '72)} can be freely multiplied to the correlation kernel 
without changing the correlation functions. Taking F{£_',ri') = i^rj ja — irj'^ /ia^ 
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and calculating 



and 



/<t(wo) 



2cr2 3cr3 



'a 2a'^ 3(t3 



C2V2 



o 



o 



Equation p.7p gives the equivalent kernel 

i(?l+«2)-i(»)?+r,i) 



e 2 ^ 



ct3/2 



(l0gfT)3/4 
rr3/2 



0F(l + f(T3/2(61oga)i/4(^,_^2)) 
A'/p2(Ci>C2) as fJ ^ cx). 



0(1) 



It is easy to verify that M^(C2,Ci) = Af<,(C2,Ci) for all ((1,(2) e K^. To prove 
weak convergence and convergence of the last particle distributions it suffices by 
Lemma [2?T] to show that there is a function Boo, which is integrable on {£,0,00) x M 
for every ^0 G such that i?oo(C) ^ \Ma{Ci 01 for all sufficiently large a. Now 



47^5/2 



-i((T«-J))^ + |(t>^+u^)+iJ(«+«)-i(cru+))i)^ 



e 2 



< 



47r5/2 
exp {(52^7 



e 2 



(cr^+25)(t^ + s=)+cri7(t-s) 



dudv 



2 I 2 c3 2(5?)^ 
+ 3" ~ 



\s + t + 25i\ 
26^} 



-dtds 



47r3/25(cr2 + 2(5) 



so choosing (5 = l/2aa gives 
M.(C,C) < 



Jt2„2 I 2 r3 



47r3/25((72 + 2^) 

exp I ^ — 



exp 



< 



(3.9) 



Sco(e,r?) 



V6 log (T 

— ^3 — 



exp 




for all sufficiently large a, which concludes the proof. 



□ 
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4. Proof of Theorem 12.51 



The proof is organized as follows: In Section [4TT] an integral representation formula 
for the sum of products of Hermite polynomials appearing in l|1.5p is derived. The 
resulting representation of the correlation kernel and a general discussion of the 
saddle point arguments used to calculate its asymptotics are given in 14.21 Sections 
14.41 and 14.61 contain the proofs of part (i) and (ii) of Theorem 12.51 respectively. The 
two cases of the proof are preceded by Sections FOl and [431 respectively, with a series 
of lemmas estimating the integrals appearing in the correlation kernel for the two 
cases. 

In the asymptotic estimates we use the following notation: C and Ci, i = 1, 2, . . . 
are generic positive constants, and the occurrence of the same symbol in different 
chains of inequalities need not denote the same number. If > and a;„/y„ 
as n ^ oo, we may write Xn = oijjn) or, provided Xn > 0, a;„ ^ y„. The relation 
Xn = 0{yn) means that there is a positive constant C such that < C?/„ for every 
n. If there is a positive constant C such that XnjC < y„ < Cxn for every n, we 
write X y„. In the proofs of lemmas and theorems where ( — (Ci, C2) is expHcitly 
stated to be fixed, constants may depend on 

4.1. Integral formula for the Hermite polynomials. The proof of Theorem 
12.51 essentially reduces to calculating asymptotics of the (appropriately rescaled) 
correlation kernel K^" . The key to this analysis is to find a suitable representation of 
the Hermite polynomials; because of the scale factor t*^ appearing in l|1.5p we cannot 
use the Christoffel-Darboux formula to simplify the sum, which is the standard 
technique. The formula derived in Lemma [4!T] differs from the one used in |5] for the 
special case (t„ — cn^^/'^, and it enables us to carry out a more complete analysis. 

Lemma 4.1. Let ri,r2 and r he positive real numbers such that ri < Tr2 and define 
the contours 



in the complex plane. Then, for any positive integer n and all complex numbers z\ 
and Z2, the identity 



(4.1) 




and 



(4.2) 



M 9 t r,.^ [t) = r2 + iteC 



(4.3) ^T''hk{zi)hk 




,w'^ — 2z2W2+2z-iW\—w\ 



dw2Awi 



k=0 



Wl — TW2 



holds. 



Proof. Recall that the orthonormal Hermite polynomials can be written 



(4.4) hk{z)^^^/^V¥k\Hk{z), 
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where 

(4.5) 77,(z) = (-l)V^^^gli, 

and that the renormahzed polynomials satisfy the generating function relation 



(4.6) g2«,2-«, 

for all complex numbers z and w. Note that for any choice of real r2, 

(4.7) e""' = / e'"'-^''"dw. 

Using Equations (|4.5p and (|4.7p gives a representation 



(4.8) Hk{z) = 



Lr 



On the other hand, Equation (|4.6p and the residue theorem yield, for any choice of 

n > 0, 

n^ 1 i e^-^ i/fc(z) 



27ri /I, w'^+i kl 

" 7T-1 

Combining the representations (|4.8p and l|4.9p gives, for any r > such that |rr2| > 



Ky — — 



e 



e 2 

27r^ Z,,^^ Vr^^ - wi 



27r^ Jy^^ 7r^2 V^^i/ wi — TW2 

as claimed, where the last equality follows since 



dw2dw 



dwi ~ 0, 

by Cauchy's theorem. □ 
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4.2. The saddle point argument. Identifying Cj = S,j + ir]j G C with (Cj:'7j) G 
and using Lemma |4?T| Equation (|1.5p can be written 



(4.10) i^^-(Ci,C2) = — T-^==exp 



C2 , , 



27r5/2y(i 



2 \ T 1+r 1 -r 



X (I) / dw2du;i, 

for appropriate choices of ri and r2, whenever < t < 1. (For the Ginibre ensemble, 
r = 0, the correlation kernel has the form 

and a representation 

^ fc! 27ri "^"(1 - w) 

analogous to l|4.3p gives a (simpler) saddle point argument, paralleling the t > 
case. The details will be skipped.) 

The correlation kernel of the rescaled point process Z^" is given by 

(4.11) i^;"(Ci,C2) 

where we will see that the scaling parameters 

(4.12) {an,bn,Cn) := \ — (a„,6„,c„), 

V Tn 

should be chosen as in the statement of the theorem. 
Define the analytic function 

(4.13) fn{w) = n\0gW + w'^ - CnW, 

where we choose the principal branch of the logarithm. We expect c„ to be close 
to (1 + Tn), the rightmost edge of the spectrum on the global scale, so the main 
contribution to the exponent in the integral in (|4.10p should be /n(w2) — fn{wi). 
The idea of the proof is now to calculate the large n asymptotics of K^" by a saddle 
point argument for /„. Define the shift parameter 

(4.14) 5„ :=c„- W — (1 + T„). 

V Tn 
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Provided < Xn — o(l), where 

(4.15) Xn = —i=7T^ ^2' 

which will turn out always to be the case, /„ has two distinct, real, positive saddle 
points, w+ > solving the equation 

(4.16) /;(w) = - + 2u;-c„ = 0. 
Explicitly, 

(4.17) w+ = 
and 

(4.18) w- = 

where the asymptotics are given by an expansion of the square root, 




Cl (l + r„)lv (l + Tn)"^" (l+r„)2 



(4.19) ^r--^ 
To analyze the behaviour of /„ close to the saddle points we observe that 

(4.20) C{w+) = 2 - = 2(1 - r„) + O { _ > 0, 



(4.21) /:(^^) = --(1 -r.n) + o{ _ J < 0, 

TV?, 



V^(l - Tn) 



and 



fi^Hw) = {-lf~^{k-l)\nw-^, for /fc > 3. 
For sequences (t„)^2 such that [/"(w-)] and |/"(w+)| become very small, we cannot 
ignore the third order terms in the Taylor expansions of /„ at the saddle points; more 
specifically this happens whenever 



lim 



>o, 



and similarly at w- , that is, whenever (Jn =0(1). This results in completely different 
asymptotic behaviour of the integral in l|4.10p depending on whether or not cr„ = 
0(1), and the choices of contours and parameters in the general parameterization of 
the integral, given in Lemma l4?2l will differ in the two cases. 

In principle, in view of l|4.20p and l|4.2ip . we would like to choose ri = w- and 
r2 = w+ in (|4.10p in order to pick up a Gaussian integral at each of the saddle points. 
However, since 

(4.22) rnW+-W^=j^^(l+0{Xn)). 

(1 - Tn) 
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the integral in (|4.10p will not converge for this choice of ri and r2 unless i5„ > 0, so 
if Sn — a slight modification will be necessary. 

For ri and r2 (depending on n) to be specified, parameterize the contours of 
integration 



(4.23) 



: s h-^ rie*''"Ms| <tt/0 
: t r2 + iant, t G M, 



where On and Q!„ are positive parameters. Moving out the main contributing factor 
from the integral, we get the following representation of the correlation kernel: 

Lemma 4.2. Let ri and r2 he any positive real numbers such that T„r2 — ri > 0. 
For i, w € R, define 

^rlit) :=/n(^2 + iant) - fn{r2) ~ »a«i(an6 - iKm) 

(4.24) =71 log ( 1 + — ^ j - alt'^ - iant (c„ + an^2 - ibnV2 - Zra) 
and 

<^>il{v)+i'$il{v) := -/„(rie™) + /„(ri) + ri(e» - l)(a„ei +iMi) 
=r^(l — 2cos2w) + ric'j(cosw — 1) — ribnT]i sinv 

(4.25) + i {{—n + l)v — r\ sin 2u + r\d^ sinu — ri&„?7i(cosw — 1)) 

where = c„ + o„^i and <i>^J and 5*^^ are i/ie real and imaginary parts of the right 
hand side of {4-25^ respectively. 

The correlation kernel K^" of the rescaled determinantal point process ZJ^'' can 
then he expressed 

jyr /> ^ X -ri6'„a„a„6„T^+^ f (c„ + a„6)^ - Tn(c„ + an^if 
t^n (Ci , C2) = , r,^ exp 



47r5/2^(l_r2) "i 4(1 + T„) 

bl{Tnr]\ + + an6) ' 



4(1 - T„) 2 
X exp{/„(r2) - fn{r\) + a„(ri^i - r2^2) + ibn{r\r\x + r2772)} 

(4.26) X //" F4v4=^(s,t)dMs, 

i J{|s|<7r/e}xR 

wftere 

U 07) pC.Grs n - ^^v{'^'^l((^ns) + ^^</,{ens) + E</M 

The technical part of the proof of Theorem [23] has now been reduced to estimating 
the integral 

(4.28) G(Ci,C2):=// F,^\'4^(s,t)dMs. 

J J{\s\<Tr/9}xR 
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In order to estimate the s-integral, we need to establish that the real function 

*J'ri(^) = ''ill — 2cos2u) + ric'j^{cosv — 1) — ribnrji sinv 

decreases monotonically as one moves away from its global maximum. This is the 
content of the following lemma. 

Lemma 4.3. Let ^ e be fixed. Put vq = —bn'ri/{c'n — 4ri) and choose ri = w_ if 
cr„ — > oo, and ri = W-{1 — e„n~^/'^) otherwise, where n~^^^ <C en = o(l). Then for 
every sufficiently large n there is a p„ with 

f Pn < ^17^7173 if'^n = 0(1) 

[ Pn < C2:^y=T;J^ otherwise, 

such that if vi £ (— 7r/2, tt/2) \ {vq — Pn, va + Pn), V2 G ti", tt] and 

\V2 - Vq\ > \vi - Vol ; 

then 

> <(«2). 

Proof, ^rii"^) is differentiable everywhere, so any local extreme points v satisfy 
^L'iv) = 0, or 



2 



COS V 

ri. Now 



(4.29) sinw ( 4ri 

provided rj 0. Put x 

(4.30) r,{c'„-4r,) 

= n(l - T„) + ((5„ + 2x{3 - 1/t„)) 
if n is sufficiently large, so 



bnV = 0, 



2 



Ox' 



O (a„V^) > 



47-1 



cosw 



> 4 - 4ri = O (a;) if cr„ 



0, 



and 



In both cases, 



4ri 



> c'„ — 4ri = O (Vn(l ~ Tn)) otherwise. 
0(1). 



4ri — c^/ cosi; 

Therefore, for any solution v to (|4.29p . sini> = O (wq) = o(l). Consider first the case 
that \v\ = 0(1). Then there are numbers ci and €2 with |ei| < jw'^l/G and |e2| < 
such that l|4.29p can be written 

{v - ei)(4ri - c'„{l + 62)) - 6„?7 = 0, 

which implies that there is a constant C such that 

W e (vO - Pn,VQ + Pn), 
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where 



Similarly, if |u — 7r| = o(l) a Taylor expansion of l|4.29p gives the condition v g 
(wo — Pn,vo + Pn) where vq = tt — bn'r]/{c'^ + 4ri). It follows that strictly 
decreasing on [wo+Pn, "Co— Pn] and strictly increasing on [vo — 2'!r+pn, VQ — pn]. Clearly, 
any possible local maxima of ^ij^{v) near vq are smaller than, say, $^^(±7r/2), so 
evaluating the order of magnitude of p„ and remembering that <i>^^ is 27r-periodic 
gives the conclusion. □ 

4.3. Estimates on G(Ci,C2) when (t„ tends to infinity. The two cases, depend- 
ing on whether or not f7„ = 0(1), require different choices of the parameters. In this 
section, estimates on 0(^1X2) are provided for the case tT„ 00. In this regime, 
the contribution to the exponent in F^l'r2 from the smaller terms depending on the 
77- variables is not negligible, so the main contribution to G(Ci,C2) actually comes 
from intervals /„, J„ which in general do not contain the saddle points of /„. 

First, the main contribution to the integral is estimated in Lemma FOl then it 
is shown in Lemma l475l that the remaining contribution is negligible, and finally in 
Lemma [131 the main contribution is explicitly evaluated. These estimates provide 
the basis for the proof of part (i) of Theorem 12.51 in Section 14.41 

Throughout this section we will suppose that f7„ ^ cx) as n — > 00 and that 



(4.31) Sn X V(l - T„)logcr„ ^ n 2/^cr„Vlogc^n> 

/(l-r„) _ 



V log cr,i VlogOn ' 

n(l-r„)^\'/'_ n-y^al^' 



(4.32) 
and 

(4.33) 6n^, , , ,1/4. 

V logCT„ J (log cr„) ^ 

Since Sn > 0, we may choose ri — w- and r2 — w+ by l|4.22p . Now choose the 
parameters an and 6n so that the exponent in becomes of order — s^/2 + 

i^/2, namely 



and 



.3/2^ 



(4.35) On = = = ^ 1 - "" + O {xf 

y'2w-{w+ ~ W-) ^n(l-T„)\ V2n(l-r„)2 



We observe that by definition the parameters are related through the equation 
(4.36) al = 1/2 + dlwl 



EDGE SCALING LIMITS FOR NON-HERMITIAN RANDOM MATRICES 



21 



and that 
(4.37) 



3/2 e 



First we estimate the contribution to G'(Ci,C2) from near the (sUghtly shifted) 
saddle points. 

Lemma 4.4. Put So = — 6'„w_6„77i and to = — «n^n??2 and let M„, such that 1 ^ 
M„ ^ \/log (T„ , 6e given. Define the sets In = {s : |s — sol ^ Afn} a^ii '^n = {t ■ 
{t-tol < M„} and let 0„ = - w_)"3/24. T/ien, /or every fixed {(1X2) e M"*, 



e 2' 



-0(1) 



-dMs 



Proof. We expand the exponent in F^^2'^^^{s,t) and estimate for s & In and t G J„. 
Note that, since 



W-{c'n~W^)el = 1 + 



a„Ci 



(c„ - 4ii;_ ) 
the expansion of the exponent becomes 

^t{t) + <^i\{ens) + t^i\iens) 



= 1 + 



1 



. ria: 



n +3 



3w3 



r + ^- 



(4.38) 



6 

6! 



2 

(1 - 
log cr„ 



."(l-T«)^ 

The imaginary second and third order terms in l|4.38p are in general not small in 
In 'x Jn, however they are almost constant. For the term involving r]2, 



(4.39) 
Since 

and 





( inalt^\ 


i(pnV2 - 


{ 3wl J 



< -^MnUY&yi{tl,Ml) < C 



V^(l-r„)5/2- 



w_(8w_ - 4) = n(3r„ - 1) + O 



^^^I^ = 0„(l + O(2;„)), 
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the terms involving r]i can be similarly estimated; 



(4.40) 



2 6 



- in^^JiK'i + (./.,)' (.(3.. + 0(1))) ^ 



< c 



V^(l-T„)5/2- 



Furthermore, 

(4.41) w-e'*'"'' - Tn[w+ + iant) 

= + i9„W-S — Tn{w+ + iant)) ( 1 + O 



y/n{l - Tn)Sn 



Using Equations l|4.38p . (|4.39p . (|4.40p and (|4.4ip and noting that the error terms 
they give rise to are all small, gives an estimate 



-i0„(ii2-»;i) 



// F<-:%^{s,t)dtds 



6 2'- 2- 



<\Rn 



K7_ — TnW+ + iW-9nS — ITnant 

2_J_ 

6 2" 2 



-dMs 



-t^ — ^s^ — iant{an^2~'ibnri2)-\-iw—9ns{an^i -\-ibnrii ) 



= (l + o(l))|i?„| 



W_ — T„Zi;+ + iW-OnS — ITnCtnt 



dMs 



|w_ - T„W+ + iW-9nSo - iTnantol ' 

where i?„ — o{l) and the conclusion follows. 



□ 



Next, it is shown that the contributions from outside the set /„ x J„ to G(Ci, C2) 
and its approximation, given in the previous lemma, are negligible. 

Lemma 4.5. Define sq, Iq, M„, /„, and J„ as in Lemma \4^ and let 0~'^T denote 
the set (-7r/6'„,7r/6l„). Then, for (Ci,C2) G K"* fixed, 



(4.42) 



(e;r^TxR)\(/„x j„) 



Fi':%{s,t)dtAs 



0(1) 
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^—■^t'^ — ^s'^ — iant{an^2—^f>n^2)+i'W-dns{an^l+ibnr]i) 
-0(1) 



dtds 



g 2 V^O^^O^ 



\w- - T„w+ + iw^0nSo - iT„a„to\ ' 



Proof. We prove the estimate l|4.42p ; the second assertion follows by the same argu- 
ment. Since 

(^-^T X M) \ (/„ X J„) = Ai U ^2 U A3, 

where = /„ x (M\ J„), A2 = ii0-^T) \ /„) x J„ and A3 = {{0^,^^) \ In) x (M\ J„), 
we consider first the integral 



(4.43) 





/ 




Jr\.j 




f 




K\J„ 


<e<;/2 



At 



cxp |n log 


1 1 iOint 




|w- + 


'.w^e^s - r„(w+ + iant)\ 



dt 



\u\>M„ \w- + iw^9„s - Tn{w+ + ianiu + to))\ 



rdu, 



in order to estimate the integral of i^^^j^i^,^ over Ai. Fix s e M and suppose first 
that »72 7^ and that 

h 

\iW-6nS - ITnanto] 



(1 - rn) ' 



this will generically be the case. Note that 



» TnW+ — W- 



(1 - Tn) 



> a„M„. 



(1 - rn) 

We can therefore choose Nn such that 

q;„M„ < a„7V„ < |w_ antol < a„iV^, 

for example 

N = ^ 

" (log a„) 1/6 
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will do. To see that the contribution to Ijg{s) from |u| > Nn is negligible, note that 



^/e„>|«|>Af„ 1^- + iw-OnS - T„(w+ + ian{u + to))\ 



du 



< 



{w+ + ian{u + to))\ 
suPAr„>|„|>M„ + iw^OnS - T„(u;+ + m„(u + to))|} 



N„<\u\<L, 



TnW+ - W_ 

e-"'/2dM 



M„<|«|<Af„ 



Equation (|4.43p therefore becomes 



p-Ml/2 



du 



*o/2 



<c- 



e 



Mn \w- + iW-9nS - Tn{'W+ + ianto)\ 



The (simpler) cases when 772 = or |iw_0„s — iTnanto] <€. bn{l — Tn)~^ can be 
handled similarly, and lead to the same estimate. Essentially the same argument as 
in the proof of Lemma [T4l now provides an estimate 



(4.44) 



F?iiMs,t)dtds 



< C- 



Mr,. ' 



Ij-(s)e 'ds 



^ |w_ - T„W+ + iw-9nS - iT„anto\ 



ds 



< C 



Now we turn to the integral over ^2. Consider, analogously with the previous 
estimate, the case that ?7i ^ and 

\iw-.9nSo - iTnant\ > T„u;+ - w^. 

Choose Nn as before and let Ln — {n{l — t„))^/^, so that L„ ^ Nn, but 6'„L„ <C 1. 
Then for t fixed, we proceed to show that the main contribution to 



(4.45) 



ds 
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comes from the intervals Af„ < |u| < iV„. First, by Lemma l43t 

r *i,i_(e„(«+so)) 
7l„<|„|<^ |w„e*e„(«+so) _ rn{w+ + ia^t)\ " 

/ n g*^i_(en(n+^o)) \ 

^supi^>l„l>M„{k- + iw-On{u + So) - T„(u;+ + ia„t)|}e*"'-*-'^"-' 
/Af„<|«|<L„e""'/'dw 
dw 

„<|u|<^ TnWjf. W— W— C0S(6'n (w + So)) 

<Ci^7^"''^°'''" < - r„) loga„)^/%.e-^"/3 « 1 



But for |s| < L„, we also have an estimate 

(e„(M+so)) 



^■i>;ii_(e„(«+so)) 

rdlt 



-1 



< 



suPm„<|«|<7V„ + iw^e.n{u + So) - T„(w+ + ia„t)|} 

JA/„<|m|<W„ 



\/Io^ 

Consequently, (|4.45p becomes 



<Ci / ^ ;^ --ds 



<C 
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SO 



(4.46) 



11^ Ft:^l{s,t)dtds 



So/2 



H^a (t) 



\W- - TnW+ + iW-dnSo - ITnttntl 



dt 



< C 



Clearly, 
(4.47) 



Fi':^l{.s,t)dtd.s 



0(1)- 



Mn - TnW+ + iW-OnSQ - «T„ari,to| ' 



which together with (|4.44p and l|4.46p completes the proof. 



□ 



We close this section with an explicit evaluation of the main contribution to the 
approximating integral from Lemma 14.41 

Lemma 4.6. For any fixed (Ci, C2) G K'*, 

exp {-\t^ ~ - iant{an£,2 - ihnm) + iw-9ns{anii + i6„?7i)} 



?«_ — T-nWj^ + iW-dnS — ITnOlnt 

27rexp {-\a\(an£,2 - ibnmY - 5^?jW^(an6 + ibnTjiY] 



dtds 



(1 + 0(1)). 



Proof. Since 5„ > 0, we have r„w+ — w_ > and the denominator can be written 
1 



W_ — TnW^ + iW-9nS — iTnUnt 



/•oo 

/ exp{-r(T„w+ - 
Jo 



W- — iw^OnS + iT„a„t}dr, 



so, using Fubini's Theorem, 



W_ — TnW^ + iW-6nS — ITnant 

(flnCl + ibn1]2f 



-27rexp<^ — ^(-a„f2 + iOn»y2) ^ — 



/•OO 

/ exp {-Air^ - (B„ + iCOr-} dr, 
Jo 
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where 



A„ = rial + wlel ^ + + o(l)), 



(1 - Tn) 



(l + o(l)), and 



Putting e = i?y^j4„^^^ and integrating by parts gives an estimate 



/•oo 

/ exp {-Anr^ - {Bn + iCn)r} dr 
Jo 



1 



B„ + iCn 



< 



< 



Bn + iC, 

An. 



■ / 2r exp {-Anr^ - (Bn + iCn)r} dr 
Jo 



Bn + iCn 
1 



2re-^^''dr 





2eAr, 



2re-^"'^'dr 



+ e 



(4.49) 



< C 



Bn + iCn\ \ Bn 



\Bn + iCn 

which, inserted into l|4.48p . gives the conclusion. 



□ 



4.4. Proof of Theoremimi part (i). To be able to apply LemmaOto (Z^")^^!, 
we first calculate the point-wise limit of (a kernel equivalent to) i^^", using the 
estimates of Section [473} and then prove that, for any real ^o, there is an integrable 
function on {£_q,oo) x M which dominates the functions {Kn"{CX)}'?^=i- 



Proof of Theorem \2.5\ (i). Suppose cr„ — > oo and choose q;„ and Qn as in the previous 
section, that is. 



and 



2(w+ — w_) 
1 



The choice of scaling parameters is given by requiring that the rescaled kernel K^" 
have a non-trivial finite limit. To be able to control the error terms we assume from 
the outset that the orders of magnitude are correct, namely that the conditions 



(4.50) 



5n ^ y/il - Tn)loga„ 



(4.51) 



log CT„ 
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and 

'n{l - Tnf 



(4.52) 



log CT„ 



are satisfied. Choose Mn such that 1 ^ M„ <C VTogCT^. It then follows from the 
estimates of Lemmas 14.41 14.51 and 14.61 that for any fixed (Ci)C2) G I^**, 

_ 27rexp {-\ai{an£.2 - ih^mf ~ \Slw1{a„^i + 

r„w;+ - + ra^ (a„^2 - «fenf?2) + 91w'^„{a.n£,i + ?6„?7i) 

Since r„?i;4- — w- > whenever (5„ > 0, we may choose ri = w_ and r2 = w+ in the 
representation (|4.26p derived for the correlation kernel of Z,^" . Recalling (|4.36p and 
(|4.37p this yields, after some simplification, 

(4.53) i^,:"(Ci,C2) 

anbnTnW-Onan exp |rilogT„ + "^lli^^^) + fn{w+) - /„ (w_ ) | 

27r3/2y'(l - tI) {t„w+ -w^+ ra^ (a„^2 - «fe«??2) + 6'2i(;2 (q^^^ + ibniji)) 
X exp |- (^(^1 + 6)(1 + o(l)) - 23/2jl(?!\)3 (^i + '^2)(1 + 

-P'ri(C2) ( 2 f 1 "riA />2 , ^2' 



where -FnlC) = exp{— i6„?i;_?7 + ibnO^upL^r] + i(j>n7f} and the small o terms in the 
exponent are both O ((5„n~^/^(l — r„)~^). The correlation kernel of a determinan- 
tal point process is not uniquely defined; clearly all correlation functions remain 
unchanged if K^"((^i, C.2) is replaced by 

(4.54) kUCi,C2) = ^^,K:HCuC2). 

In order for the exponent in (|4.53p to have a non-trivial finite limit depending on 
the variables Cj, we choose a„ and 6„ such that 

C4 t^^l lim - lim ^"^ ^"'^'^ - i 

^ ™ (1 - r2) " ™ 23/2^^(1 - T,03 " 2 

as rt ^ 00. To evaluate the constant exponential factor in l|4.53p the precise asymp- 
totics (|4.19p are needed. Expansions of the logarithmic terms, and a considerable 
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amount of subsequent algebraic manipulation, give 

4(1 + T„j 

4(1 - Tn) ,22 / N , \/2T„n(5„ 

j;2 ( '''ni^'i — 2t„ — 2t^ — 2t^ + t!^^) \ , q f 



4(l + r„){l-T„)3 ; VV^(l-^n)5 



After an asymptotic expansion of the denominator of the constant factor in l|4.53p , 
Equation (|4.54p thus becomes 

(4.56) <(Ci,C2) 



27r3/2y(rT^^„ (l + zi^(,;i - r;^)) 

The exact choice of (5„ is now given by requiring that if„(Ci, C2) have a finite limit, 

say 



™ 27r3/V(l + ^.)X; ~ 7^' 

or equivalently, in view of l|4.55p . 



(4.57) lim + nV^(l - r„)^^-^/^ exp ( ^ = 1. 

Equation (|4.57p is satisfied for 



" ^ ^ ^^It'''^^ " Tn)61ogg„ 



so by l|4.55p we may choose 



' {1 + Tn) / (l-r„) 
T„ y 61og(T„ 



and 

V 6l0g Cr„ 



for! 
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which are the choices in the statement of the theorem. Since these parameters satisfy 
the assumptions (|4.50p through (|4.52p and ^ 

e-^(Ci+?2)-i(')?+r,2') 

i^JCl,C2)- 7 ^(l+0(l))^Mp2(Cl,C2), 



as 71 00. 

It remains to prove that, for any given <^o G Kn"iCX) is dominated by an 
integrable function on {£,q,oo) x R for every sufficiently large n. To estimate the 
integral 



'(e;7^T)xi 

from above, let e < 1 be fixed and consider any n so large that 

(4-58 Ko < e 77- T72 n 

For such a choice of n, it follows that 

3a^a„|^o| ^ 4(1 -T„) 



W+ (1 + Tn)bl 



£2 



Cn - iW- (1 + Tn)bl ' 



and that 



(r„w+-w_) '(1 + T,i) 



log CT,, 



To capture some of the oscillatory terms in the exponent of i^^'f , we may again 
change contours of integration, by Cauchy's theorem. Replacing the real line by the 
contour E — za„^/(2a„) in the t-integral, and putting t' — t + ianS,/{2an), gives 



(4.59) 



// i^£^,^^(s,t)dMs 



< 



{t„w+ - w_ - r„a„|^o|) J(e^iT) 
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Now, using that f^{w+) = and the inequality log(l + x) < x gives an estimate 
/ exp {t' - ia„C/(2a„))| dt'l 

< ^ cxp {Re (s« ^ {t' - ia„^/(2a„))) } dt' 



< 



e 2al - 



naz 



2 V " iw++an^/2y 



+nlog ( 1 + ] + a„5(w+ - c„/2) 



+ anbnVt' 



(4.60) 



< exp I - (1 - ei) + a„6„7?t' - ^ | 



To obtain a uniform bound in n of the first integral, recall that 

<-4»_>(c„-4»_,(l-J^ 

where = c„ + a„^. Suppose without loss of generality that ry < 0. Then there 
aVri & [0,7r/2) such that (u^) = and since 

(ii,,) = Aw^_ cos2t)^ — w_c^ cosi;^ + W-bnijsinvri 
(4.61) < w_(4u;_ - c„) < 0, 

a saddle point argument gives 



(4.62) / 



e as 



< Cie 
Consider the difference 



exp 



1 / c' - 4w_ 



2 V c„ - Aw- 



s — 



ds < Ce 



2(1-62) 



2(1 - 62) 

(1 - €2) sin^ , 



> 



+ bnW-rjsinVr, — W-c'^{cosVr, — 1) — ■w?.(l — cos2t;^)) 

c^ui- (cos fr, — 1) — w'^{l — cos2t)^) =: g{vn). 



Since 



201 



g'{Vr,) = CnW- SmVr, ( — — (1 — £2) COS t;,, ) > 0, 
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g is increasing on [0, 7r/2) so, noting that g{0) = 0, it follows from (|4.62p that 

This, together with the estimates l|4.59p and (|4.60p . inserted into (|4.26p shows that 

\Kl^iC. 01 <C. exp { {^-^^^ - - ^-)) e - ^ 

2 V2(l-T„) (l-6l) (1-62)^/ 

+ 2 (l-^2)j/ 

(4.63) <Q„e-(i-^)«-(i-')''', 

where C^^ denotes a constant depending on ^o- This provides an integrable bound 
on ((, 0, so by Lemma (|2.ip and the point-wise convergence of if^ to Mp2, 
converges weakly to Zp and the last particle distribution F,'^" converges to Fq. □ 

4.5. Estimates on G(Ci,C2) when cr„ tends to a finite limit. Throughout this 
section, we will suppose that (t„ cr G [0, oo). Again, we prove a series of estimates 
on G(Ci, C2) for this case. First, in Lemma [4?7l on the main contribution from close 
to the saddle points, and then in Lemmas 14.91 and 14.101 on the remaining, small 
contributions. 

Fix the choices 5n — 0, a„ = V27v^n^^/® and 5„ = \/2t„(1 — t„) of the scaHng 
parameters. Since 5n — 0, w- — t„w+ — and we must choose contours slightly 
removed from the saddle points in order for G(Ci,C2) to converge, say ri = w-^ — 
TnOin^n and r2 = + Q;„e„, where e„ <C 1 will be specified. The quadratic terms 
in the expansions of /„ at the saddle points w_ and w+ may now be arbitrarily 
small (depending on how small cr„ becomes), so we choose a„ = (2r„)^^/^n^/^ and 
9n — n"^/'^ so that the third order terms become of order one. 

To simplify the calculations by using the fact that w+ is a saddle point of /„, we 
note that, by definition, 

(4.64) F^lj-{s,t) = i^Ci;C2^(s,t _ ie,Je-^"('^^)+/"('"+)+""^"(""«=-^''"''^). 

We begin by approximating the main contribution, from close to the saddle points. 

Lemma 4.7. Given e„ and T„, choose ri = w- — T„a„e„ and r2 = w+ + and 
define the contour 

In ■ (-r„,T„) — > c, 

ln{t) = t + ien- 

Let {(1X2) G be fixed. 

(i) If a > or an = 0{n^^^^^) choose Tn — n'^^ and e„ — n^'^^ for some 
< fci < 1/15 and h < k2 < 2fci. 
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(ii) If a = but cr„n^/-^^ oo , choose Tn = cr^"^^ and e„ = cr™^, for some 
< mi < 1/2, and mi < m2 < 2mi. 



For these choices. 



(4.65) J^ni/6e-^(''?+''2) ff F^liUs,t - ien)dtds 

/ 77 ^ dudv + o(l). 



Remark 4.8. The case in which a = but condition (ii) is not satisfied causes no 
problem, but for simplicity we omit the details. 

Proof Note that in both cases, 1 < r„ «; n^^^^ and < e„ < T~^. Put 
i =t — icn and s = s + icn- Using the expansion 

4>^i(w) + i*^i(t;) =v{-ribnm +i{l-n- 2r? + nc'J) 
- 2(^i''« "^rjw -2 

6 6 

with ri = w_ — T„a„e„, and keeping in mind that e„ ;:?> n~^/^^, one obtains, after 
some messy but straight-forward calculations. 



<c(e„ + n-i/V). 
Using that f^(w+) = 0, the expansion of the exponent in F^^^^+i^'^ ~ ^^n) becomes 
(4.66) ^% (i) + {0ns) + z^il (0„s) 

Note also that 



ii(-T„,Tn)M*(S ~ 'J ~ ^£"1 
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Consider first the case that either cr > 0, or cr„ = 0(n ^Z^^). Then, by (|4.26p and 



e 2 



i(cr„t+r;2)^-|t''-4iC2-5(CT„s+j)i)^ + is^+is^i 



g-i(cr„t+Jj2)^-|P-it{2-|(o-„S+»)i)^ + 4s^+i5?i 



-dMs 



<C: 



2 n- 



-l/3rp4 



i{s — t) — 2e^ 



/(_T„,T„)2 - - 2e„| 

<C3iog — fe-i/3 + r„e„) 

(4.67) <C\ogn{n^^^~'^/^ +n^^-^^^ 
Next, suppose n^"^/^^ < ct„ = o(l). Then 



"^1/6 



ri ,10+ 

(-T„,T„)2 
g— f f — it{2 + f +isii 

I 7 ^ o dtds 

'(-T„,T„)2 «(s — — 2e„ 

<Ci (a„T„ + n-i/3j^4^e„) 



i{s — t) — 2e„ 



dMs 



(4.68) 



<C3 (a„r2 + T^n-i/s + r„e„^ log 

<ClogCr„ (^(T^ 



(-T„,T„)2 \i{s -t)- 2e„| 



dMs 



2mi ^ 5mi^-l/3 _^ cr-'"i+'"^ 



Letting u = s + isn and i; = — (t — ie„), l|4.67p and (|4.68p give the conclusion. □ 

We turn to an estimate of the contribution to the integral of F^^'^^isjt — if-n) 
from outside the set (— T„,T„)^. 
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Lemma 4.9. Let Tn and e„ be as in Lemma \4-7\ and put ri = W- — TnCXnen- Let 
{C1X2) e 1R4 be fixed. Then 



(4.69) // 
Proof. Note that 



((e;:^T)xR)\(-T„,T„)2 



exp {SC?^ (t) + (0„s) + iM/Cj (0„s)} 



(4.70) 



{w- - T„Q!„e„)e*^"* - r„(w+ + za„t) 
<^_gRc(H<2(f)) cE,ci(e„s) 



so we estimate the integrals over s and t separately. 

Fori = 0(ni/3), we can estimate the logarithm in (i) from above by its third 
order Taylor polynomial, giving 

n . 



Re 



(Ei^^ (i)) =1 log (1 + n-y^el + t^) + 2.„n-i/3) 



e„ri 



„i/3 

2/3 _ . /+2 ,2 



2t„ 



< - 



1 / (T„< 



(4.71) 
where 



^-1/3^4 



T„r/2 



3„2/3 



1 2 .2 t ''^n^ri ^ 

2 ^T„ 



t2 t2 



To control Qn{t), we observe that if e„n^/3r„ ^ > 1/4, we can write 



Qn{t) < 



1-Ci- 



i2/3 t6 



C'2e^, 



which is either negative or of order o(l). On the other hand, if tnn^^^T^ ^ > 1/4 we 
have the estimate 

2t^ 



At)< 



„-l/3^4 



1 



3n2/3 



(1 + 1/3) < 



In both cases it follows from l|4.7ip that 

Re(s£^^(t)) 

for \t\ < v}/^ , so 

Ro(h^2 (t) 



(4.72) 



+ ''')dt < C- 



T„<|t|<ni/3 
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We turn to an estimate for \t\ > n^^^. Changing variables u = n ^^^t and letting 
5 = {e ~ 2)/(2e) > gives an estimate 



J\t\ 



'|t|>ni/3 



=ni/3 / exp {| (log (1 + ^.2+ „-2/3,2 2e„n-i/3)) 



|«|>i 

anU ' UT]2 - ' + 6^ - tn^^ ) clw 



j2 4/3^2 



<2ni/3 ^ exp || (^log(l + u^) - (1 - <5)u2 + ""''^ 



i2/3, 



2 



<Cini/3 / exp f 5 (log(2M2) - (l - 5)u^) } du 
■ 'u>i ^ 2 J 



2 J \n{l - 5) 

6 + e 



^4 + 2e 

(4.73) <C?i-i/^e"''", 

for some fc > 0. Similarly, noting that <I>^i(0„s) < -e„s^ + Ci (n^^/^s** + l) and 
recalling that e„ > n"^/^^, 



f e^'l^'-'^dsKcJ 



'Tr,<\s\<-^ JT„<\s\<^ 



"'T„<|s|<ni/12 



■,,1/6 



e 

T„<|s|<ni/12 



(4.74) <cLv3,-"^'™ + £r;^y 

It is also clear from (|4.72p and (|4.73p that 

^K-^(*))dtx /'e-"*^dt<Ce-V2 
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and similarly for the s-integral, so by (|4.70p . 
which concludes the proof. 



((e^^T)xB)\(-T„,T„) 



5/2^ 



□ 



Finally, we show that the contribution to the limiting integral from outside {—T„ , TnY 
is small. 

Lemma 4.10. Let e„ and T„ he as in Lemma {4-'^ - Then, for any fixed (Ci,C2) € 
ib4 



Proof. For any u — s + ie, we can estimate 

'7\7n 



-dudv 



o(l) 



dv 



< 



i{u + v) 

|i|>T„ z(s + ien + t + 

exp |-i((Tt - 772)^ - e„t2 + + ~ ^2en| 



dt 



<Ci- 



\t\>T„ 
1 



e"^"* dt < C- 



|t|>T,. 



2e. 



f2T 



dt 



Therefore 



7\7n Jl 



e 2 



z(w + v) 



-dudw 



< C 



5/2„ 



and by the symmetry of the variables u and v the conclusion follows. 



□ 



4.6. Proof of Theorem l2.5L part (ii). Using the estimates of Section |43] to prove 
the point-wise convergence of (a kernel equivalent to) to , and then proving 
that -ftr^"(C,C) is dominated by an integrable function, we can apply Lemma ETT] in 
this case as well. 

Proof of Theorem \2.5\ (ii). Suppose cr„ ^ ct G [0,oo) as n 00. Let (5„ = 0, 
a„ — (2r„)^^/^n^/^ and 9n = rT^I'^ . In order to tidy up the calculations slightly. 
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we will actually prove the theorem for the choice of parameters a = T„n~^/^ and 
b — T„n^^/'^cr„, differing by a factor r„ from the choice in the statement of the 
theorem; this clearly makes no difference in the limit. Thus a„ — y/2T^n~^^^ and 
bn = •\/2t„(1 — Tn) as defined by i|4.12p . Given e„ and r„ as in Lemmas 14.71 and 1 4.91 
choose ri — w_ — r„Q;„e„ and r2 = w+ + a„e„. Let (Ci, C2) G I^'* be fixed. 
By the estimates of Lemmas 14. 7( 14.91 and 14.101 it follows that 

J J{-Tr/9,Tr/e)xR 

= -M^„(Ci,C2) + o(l). 
Inserting into l|4.26p and comparing with the kernel Mo-„ , results in the expression 

^;"(Ci,C2)/Af.„(Ci,C2) 

= exp|ir„Vn(l -r„)(?7i -773) + 2("i ^ ) (^2 - ^-n^i) 
+ ^((1 " r'jr^l + (1 - + l^^V^ | (1 + o(l)) 

= exp|iT„^n(l - Tn){r]i - ?72)| (1 + o(l)). 
Considering the equivalent kernel 



Ji^^"'(Cl,C2) = exp{-iTny/n{l - T„)(?7i - 772 ) ji^^" (Cl , C2 ) , 

we see that 

i^;"'(Ci, C2) = M,„(Ci, C2)(l + 0(1)) ^ M,(Ci, C2) as n ^ 00. 

It remains to prove that Kn"{(,() is dominated, for every sufficiently large n, 
by an integrable function on {^0,00) x R. This part of the proof is similar to the 
corresponding step in the proof of the first part of Theorem 12.51 in Section 14.41 
although the subtleties of that case are not present here. 

Let e < 1 be fixed and consider any n so large that 

(4.76) l^ol < en'/\ 

Choose 



and 



7-2 = J — (1 + ei), 

V ^■''n 



where ti — {\ + a'^)n To estimate the integral 

II F^ii^{s,t)dtds, 
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replace the real line by the contour M + ja„(l — r„)^/(4Q!„) when integrating in t, 
which is allowed by the condition (|4.76p on n, and put 



t' = t + iS^ ■- t + i 



4a„ 



This gives 



(4.77) 



(e^^T)xR 

< 



(e;:^T)xR 



(V2T„nei - T„a„5|^o|) Je-^i 
Using the inequality log(l + x) < x provides an estimate 

gRe(H«^(t'-»5e))di' 

,2 f3al + 2-2{l + al)e/r2 



< / exp < —t 



- CTnVt' 



2r„(l + ei) 

-a„(S(a„ - Q;„(5)f^ + a„(5 - c„ + 2r2^ dt' 
(4.78) <C exp 1^(1 + e) - a„5(a„ - a„(5)f^ + a„(5(n/r2 - c„ + 2r2)^ 



Suppose without loss of generality that 77 < 0. Then there is a w,, g [0,7r/2) such 
that $£_'(w^) = and since 

i^v) —^Ti cos 2w,, — ricj^ cos + ri6„?7 sin 
<4r? - ri4 

< - n^/^crj; - ein(3r„ - 1) + V^Onl^ol 
<- 71^/3 (l + 3r„ (7,2 ) < 0, 

a saddle point argument shows that 



|e„s|<7r 



< Cie 



exp 



(ric;-4r?)02 



(s - w,,/6l„) Us < Ce 
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Consider the difference 



6 



1 



Since 



{jf + 6ri6„?7sinu^ — 6ric^ ((cosu^ — 1) + r\{\ — cos2w,j))) 
3 

> - 2^1.^1 sa\^{vri) - ric^(cosu^ - 1) - rl{l - cos2-d^) =: g{vr^). 



?'K) =r2(4-362)sinv„ 
>r^(4- 36^Jsint;^ 



ri(4-362) 

C,i - On 1^0 I 



4ri 



cos Vrf] > 0, 



(7 is increasing and positive on [0,7r/2), which shows that 
(4.79) / e^'i^^-^^ds^Ce'''/^ 

Recall that r2 is "almost" a saddle point of /„ , so 



2r, 



r2 



2nei((l -r„) +ei) 



< 6\/2(l + (7^)n-i 



-1/6 



Therefore, by (|4:26l) and the estimates (|4T7|) . (|478l) and (|479l) . 



(1 - T„) 1 

4 



,4(1 + T„) 16 
X exp {- ((1 + r„)(l + al) + a„(5(n/r2 - c„ + 2r2)) 
'r„(l + T„) 2 + e' 



X exp , 

^ ' ' 2 6 

which is an integrable function on (fo, oo) x R. It follows from Lemma (|2.ip and the 
point- wise convergence of K^"' to that converges weakly to Zc and that the 
last particle distribution F^"-(t) converges to Fa-{t). □ 

Acknowledgement: I would like to express my sincere gratitude to Kurt Johansson 
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Appendix A. 



For sequences of determinantal processes, the following criterion relates weak 
convergence to convergence of the correlation kernels. 

Lemma A.l. For each n > 1, let Xn he a determinantal point processes on a 
complete separable metric space A, with correlation kernel Kn- If K : ^ C is a 
function such that for every compact set A C A, 
(A.l) 

oo 

^ - / det - det C,))t=i <^'^iO ^Oasn^^, 

then K is the correlation kernel of a determinantal point process X , and A„ con- 
verges weakly to X as n oo. 

Proof. Let {A^}*^^]^ be a family of bounded disjoint Borel sets in A and t = {ti, . . . ,tk) G 
[—In 2,00)'' be given. Define the function 0t(s) = X]j=i(6~*' ~ ^)XAj{s). The sup- 
port of (p is the compact set A — u'^^^^Aj, and \<j>{s)\ < 1. The Laplace transform of 
the fc-dimensional distributions of A„ can then be written 

r r . 

^n{t): =E„ cxp <^ - ^ I A„ n Aj 



as 71 ^ 00, by l|2.2p and the hypothesis. Convergence of the Laplace transforms for all 
i in a neighbourhood of the origin implies convergence of the joint distributions. The 
existence of a determinantal point process with correlation kernel K is an immediate 
consequence of the fact that the Kn all are correlation kernels; the consistency 
conditions e.g. in [13] are easily seen to be satisfied. 

□ 

Lemma A. 2. Let Z be a determinantal point process on with a Hermitian cor- 
relation kernel K and suppose that 



i^(C,C)dC<oo 



(Co,oo)xl 
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for every G K- Then Z has a last particle almost surely and the distribution 
function F of the last particle is given by 

(A.2) F{t) = f2 ^^ / det{K{QXjm.,=id'^C 

r=0 "'((t,oo)xR)'' 

// K defines a locally trace class integral operator on -L^(R^), F(t) can be more 
compactly expressed as a Fredholm determinant 

(A.3) F{t) = det(/ - i^)L2((t,oo)xR)- 

Proof. Let t e M be given and suppose s > t. Applying (|2.2p . with 4'{s,t) = 
X((t.s)x(-s.s)), gives 

P[|Zn((<,s)x(-s,s))| = 0] 



/ det(if(C.,0))-,=id'-C. 
t^n ^! JHt.s)y.(~s.s)Y 



r=0 ^' "'((t,s)x(-s,s)) 

Because if is Hermitian, the matrix {K{QXj)Yi,j=i is positive definite, so by Hadamard's 
inequahty 



Eh/ dct(ifte,0))L=id^C<EH / 



— exp 

Since correlation functions are non-negative, 
(-1)- 



^(C,C)dC 

t,oo)xR J 
J{t,oo)xR J 



det(if(C.,0))L=id'-C 

^! 7((t,s)x(-s,s))'- 

det(K(C.,0))L=id'-C, 

SO by the dominated convergence theorem 

F{t) = lim P [\Z n ((<, s) X (-S, s)) I = 0] 

s — >oc 

= ES^/ det(if(C.,0))L=id'-C. 

r=0 "'((t,oo)xR)'' 

□ 

Proof of Lemma Wl\ We begin with the convergence of the last particle distribution. 
The existence and distribution of a last particle for Zn is given by Lemma IA.2I By 
Hadamard's inequality. 



( llldet(if„(C.,0))-,.i 



■ i=l ■ i=l 
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Now 

Ji{t,oo)xmy ^' ^' \J(,t,oo)xM ) 



and 



r=0 



SO applying the dominated convergence theorem twice gives the conclusion 



lim^ F„(<) = Jim^ I iKn{QXj))lj=i d^^C 



(-ir 



= £ / lim det (i^„(C„ d'-C 

= ^det(if(C.,0))[.,.id^C 

To verify the condition for weak convergence of point processes stated in Lemma 
lA.ll the same argument applies since an arbitrary compact set A C is contained 
in a set of the form {t, oo) x R. □ 
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